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Abstract 

This paper presents a double AR model without intercept (DARWIN model) and provides us 
a new way to study the non-stationary heteroskedastic time series. It is shown that the DARWIN 
model is always non-stationary and heteroskedastic, and its sample properties depends on the 
Lyapunov exponent. An easy-to-implement estimator is proposed for the Lyapunov exponent, 
and it is unbiased, strongly consistent and asymptotically normal. Based on this estimator, a 
powerful test is constructed for testing the stability of the model. Moreover, this paper proposes 
the quasi-maximum likelihood estimator (QMLE) for the DARWIN model, which has an explicit 
form. The strong consistency and asymptotical normality of the QMLE are established regardless 
of the sign of the Lyapunov exponent. Simulation studies are conducted to assess the performance 
of the estimation and testing and an empirical example is given for illustrating the usefulness of 
the DARWIN model. 

Key words: DAR model; DARWIN model; Geometric Brownian motion; Heteroscedasticity; 
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1 Introduction 


Time-varying volatility has been crucial in modeling economic and financial time series. Since the 
seminal work of Engle (1982), the autoregressive conditional heteroscedasticity (ARCH) model and 
its numerous variants have been widely used; see, e.g., Bollerslev et al. (1992) and Francq and 
Zakoian (2010). Among them, the first order double autoregressive (DAR) model has attracted 
much attention, which takes the form 


yt = fayt-i + m 


y4jo + ao2/?_i, 


( 1 . 1 ) 


where 4>o £ (—00,00), ojq > 0, ao > 0, {77*} is a sequence of independent and identically distributed 
(i.i.d.) random variables with zero mean and unit variance, and rjt is independent of {yj',j < £}• 
Model 0> was initially introduced by Weiss (1984), and the term ‘double autoregressive ’ was coined 
by Ling (2004) since both the conditional mean and variance functions are regressions only on the 


observed data and its square, respectively. Clearly, model (1.1) belongs to the class of ARMA- 
ARCH models in Weiss (1984) and of nonlinear AR models in Cline and Pu (2004), but it is 
different from Engle’s ARCH model if 4>o 7^ 0. Its higher-order extension and generalization can be 
found in Weiss (1984), Lu (1998), Ling (2007), Zhu and Ling (2013), Guo, Ling and Zhu (2014), Li, 
Ling and Zakoian (2015), Li, Ling and Zhang (2015) and many others. 


The stationarity conditions and the associate inferential theory of model (1.1) have been well 
studied under the compact parameter space A = {Ao = (4>o,luo, ckq) : 


< 4>,ui < coq < Co,a < 

a 0 < d}, where (j>, w, Co, a and a are some finite positive constants. Early contributions in this 
context include Guegan and Deibolt (1994) and Borkovec and Kliippelberg (2001), who derived the 


sufficient and necessary condition of weak stationarity of model (1.1). Recently, Chen, Li and Ling 
(2014) proved that when rjt is symmetric, model 0 is strictly stationary if and only if the Lyapunov 
exponent 70 = E log |0 O + rjtyfa o| < 0. 


By assuming that the true value Ao is an interior point of A, the inference of model (1.1) based 
on quasi-maximum likelihood estimation exhibits quite different phenomenon in terms of 70. For 
instance, when 70 < 0, Ling (2004) showed that the quasi-maximum likelihood estimator (QMLE) 
of Ao is consistent and asymptotically normal; and when 70 > 0, Ling and Li(2008) and Chen, 
Li and Ling (2014) demonstrated that the (unconstrained) QMLE of ((/>o,ao) is consistent and 
asymptotically normal, but the intercept term wq cannot be consistently estimated. This phenomenon 
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can also be found for the least absolute deviation estimator and robust quasi-maximum likelihood 
estimation in Chan and Peng (2005), Zhu and Ling (2013), and Li, Guo, and Li (2015). 

In all aforementioned work, the positivity of Wo is essential for the strict stationarity and (robust) 
quasi-maximum likelihood estimation, although the intercept wo is not involved directly in 70 . The 
case that uq = 0 would be meaningful but hardly touched so far. As one motivation of this paper, 
it is of interest to fill in this gap from a theoretical viewpoint. Another parallel motivation of this 


study is from application. The importance of model (1.1) and its higher-order extension has been 


well demonstrated by empirical studies; see, e.g., Ling (2004), Zhu and Ling (2013), and many others. 
In these applications, however, an interesting finding is that the intercept term ujq is often very close 
to zero; see also the illustrating example in Section 6 below. As a result, it is intuitively appealing 
to study a first-order double AR model without intercept (abbreviated to DARWIN(l)) as follows: 


yt = <t>m-i + m\! 


( 1 . 2 ) 


where all notations inherit from model (1.1) with ckq > 0. In the special case of </>o = 0, model (1.2) 


becomes the ARCH model without intercept in Hafner and Preminger (2015). In this paper, we are 


concerned with the probability structure of model ( 1 . 2 ) and the estimation and inference of 70 and 


(</>o,«o)- It is surprising that its probabilistic structures and asymptotics of the QMLE are totally 


different from those of model (1.1). This is out of our expectation. Moreover, model (1.2) is always 


non-stationary and heteroskedastic regardless of the sign of 70 , and hence it provides us a new way to 
model the non-stationary heteroskedastic time series. Specifically, it can be seen that when 70 = 0, 
the conditional volatility of yt is a nondegenerate random variable oscillating randomly between zero 
and infinity over time, while when 70 7 ^ 0, this is not the case. In this sense, model (2) is stable if 
70 = 0, and unstable otherwise; see also Hafner and Preminger (2015) for the same argument in the 
ARCH model without intercept. 

This paper is organized as follows. Section 2 considers sample path properties of yt in model 


(1.2). Section 3 proposes a new estimator of 70 and discusses its asymptotic theory and also a test 


for stability. Section 4 focuses on the QMLE of (4> o,ao) in model (1.2) and derives its asymptotic 


properties. Sections 5 and 6 investigate the numerical properties of the proposed procedures using 
both simulated and real data. The conclusions are offered in Section 7. All technical proofs are 
relegated to the Appendix. 
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2 Sample path properties 


In this section, we study sample path properties of yt in model (1.2), when ry is symmetric but not 


necessarily has mean zero and variance one. To accomplish it, we consider an auxiliary process 


xt = 


(f>0Xt-l + r) t y/o!oa^_i 


where the initial value xo = 12 / 01 • It is straightforward to see that 

xt = \(/>o + VtV®o\ x t-i or logx t = log 1 0 o + r] ty /a^\ +logx t _i. 


( 2 . 1 ) 


( 2 . 2 ) 


Note that when ry is symmetric, it is readily seen that 


{xt} = {\yt\y, 


(2.3) 


see, e.g., Borkovec (2000) and Borkovec and Kliippelberg (2001). Thus, \yt\ either converges to zero 
or diverges to infinity almost surely at an exponential rate as t —> oo, according to 70 < 0 or 70 > 0 , 


respectively. This result sharply differs from the one in model (1.1) when 70 < 0. 


Next, we precisely characterize asymptotic distribution of \yt\ after a suitable renormalization. 


Let [a] be the integral part of any real number a. The expression (2.2) implies that for any s £ [0,1], 

1 x\ ns \ 1 1 

—j= log - 77-4 - T = — > (log \<t >0 + 7?iVOo| - 7o) + —j= logXQ. 

y/n exp{[nsj 7 o} V n x V n 


By (2.3) and Donsker’s Theorem in Billingsley (1999, p.90), we have the following theorem: 


Theorem 2.1. Suppose that {ry} is a sequence of i.i.d. symmetric random variables satisfying 
< 7(7 = var(log \ (/)q + yty/aol) £ (0, 00 ). If yo is symmetric with P(yo = 0 ) = 0 , and independent of 
{ht • t > 1}, then 


\V[n 


exp{< 7 7 o lB(s)} as n —>• 00 m D[ 0 , 1 ], 


exp(s 7 oy / n) 

where '=> ’ denotes weak convergence, B(s) is a standard Brownian motion on [0,1], and D[0,1] is 
the space of functions defined on [0,1], which are right continuous and have left limits, endowed with 
the Skorokhod topology. 

Theorem 2.1 has many implications. First, it implies that yt is non-stationary and heteroskedastic 
regardless of the sign of 70 . Second, its sample path property depends on the sign of 70 . Precisely, it 
indicates that, if 70 = 0 , the conditional volatility of yt is a nondegenerate random variable oscillating 
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Figure 1: One sample path of {yt}\^\ 
respectively. 


in model ( 1 . 2 ) corresponding to 70 < 0 , 70 


0 , and 70 > 0 , 


randomly between zero and infinity over time; otherwise, \yt\ either converges to zero or diverges to 
infinity almost surely (a.s.) as t —> 00 . To illustrate this, we present one sample path of {yt} of 


model (1.2) with rj t ~ 1V(0,1), 0o = 0.5, and «o = 3.1 (i.e., 70 < 0), 3.3058 (i.e., 70 = 0), or 3.5 


(i.e., 70 > 0 ). The plots are depicted in Fig{TJ They show clear evidence of different sample path 


properties with various 70 . To distinguish them, model (1.2) is called stable if 70 = 0, and unstable 


otherwise; see Hafner and Preminger (2015) for the same argument in the ARCH model without 


intercept. Since 70 plays a key role in determining the stability of model ( 1 . 2 ), it is desirable to 
consider its estimation and inference in the next section. 


3 Estimation of Lyapunov exponent 

In this section, we propose a simple estimator for 70 . This estimator requires only the data y n 
and 2/0 but shares good properties. As in Section 2 , we assume that rj t is symmetric and not 


necessarily has mean zero and variance one. The basic idea is very intuitive. By (2.1), we have 
log(x t /x t -i) = log | 0 o + Vty/a o|, and then 


1 1 n 1 n _ 

(log log Xq) = -y^ log(x t /x t -l) = -Y] log |0O + VtV®o\ 
n n z ' n z ' 


70 a.s. 


t =1 


t =1 
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provided that 70 < 00 . Consequently, an easy-to-implement estimator of 70 is defined by 


1 


In = (log | Vn | - log 12/01) - 


n 


The following theorem states that is unbiased and asymptotically normal. 

Theorem 3.1. If the assumptions in Theorem 2.1 hold, then, (i) is unbiased; (ii) 
a.s.; and (iii) y/n{ 7 ™ — 70 ) =$■ N{0,a^ ) as n —> 00 , where cri is defined as in Theorem 2.1. 


7o 


As an application, consider a testing problem whether model (1.2) is stable or not, i.e., 


#0 : 7 o = 0 v.s. Hi : 70 / 0 . 


Based on Theorem 3.1, the proposed test statistic is given by 

T n = (3.1) 

dry 

where ^ X)™ = i{l°s(|yt|/l2/t-i |)} 2 — 7 n- Under Ho, it is not hard to prove that T n = =7- N( 0,1) as 
n —> 00 . Thus, Ho is rejected at the significance level (5 G (0,1) when \T n \ > <h _1 (/3/2), where 4>(-) 
is the cumulative distribution function of N{ 0,1). 


4 Quasi-maximum likelihood estimation 


Let 9 = (4>,a) T be the unknown parameter of model (1.2) with true value 8q = (</>o,oo) T - Denote 
the parameter space by 0 = 1Z x (0, 00 ). Assume that {y\, ■ ■ ■ , y n } are generated by model (1.2). 
When rjt ~ 1V(0,1), the log-likelihood function (ignoring a constant) can be written as 

(:yt ~ fiyt-1) 2 


1 f 

L n{fi) = 5 ^*t( 0 ) with I t (0) =--l\og(ay 2 t _i) + 
t =l ^ 


a Vt -1 


Then, the quasi-maximum likelihood estimator (QMLE) of 9q is defined as 


9n := (0n,Sn) T = argmax L n (6) . 

By setting dL n (6)/30 = 0, it is not hard to see that 0 n has a unique explicit expression with 


1 


E 


yt 


and (I n = — 


n ' yt -1 
t =1 y 


1 n 
n ^ 


t =1 


{yt - (j> n yt -1 
Vt-i 


With the help of this explicit expression, it is convenient to obtain the asymptotic properties of 9 n 
in the following theorem. 

Theorem 4.1. Suppose that Eyt = 0 and Eyf = 1. Then, 
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(i) (f) n is unbiased, but a n is asymptotically unbiased; 

(ii) 9 n — > 6q a.s. as — »• oo; 

(iii) Furthermore, if Erf = 0 and k 4 = Erf < 00, then 

Vfi{6 n -e o )=^N(0,E) 


as n —>• 00, where S = diag(ao) («4 — l) a o)- 


Remark f.l. In Theorem 4.1 (i) and (ii), the symmetric condition of r]t is not required but we 
need one more condition that Erf = 0 for part (iii) of this theorem. This is to guarantee the existence 
of covariance E. From the proof of this theorem, we have 


lim co v(rfn((j) n - (j> 0 ), rfn(a n - a 0 )) = o^ 2 Erf x 


lim — 

n —>00 77, 


1 U 

-^F^sign^t- 1 )] 


t =1 


since yt is non-stationary, the limit in the proceeding equation may not exist unless Erf = 0. 

Remark f.2. Although a n is biased, its re-scaled version a* = is unbiased. 

Remark 4-3 . Compared to the asymptotic property of the QMLE in Ling (2004) and Chen, Li 
and Ling (2014), Theorem 4.1 has two interesting features. First, the compactness of the parameter 
space 0 is not necessary. Second, the asymptotic covariances of 0 n are the same in both the stable 
and unstable cases. 

Remark 4-4- As a possible application of Theorem 4.1, one may consider a natural plug-in 
estimator of 70 as 7 * = \ Yrf=\ log \(j) n + rjtV^n\, where rf t = (yt - f>nVt- 1 )/ \Ja n y^_ 1 is the residual 
of model (1.2). However, unlike Francq and Zakoian (2012), the asymptotic normality of 7 * requires 
that E{((j) 0 + yty/orf) 1 } < CO) which does not exist generally, even for r] t being standard normal. 
See also Li, Guo and Li (2015) for more relevant discussions. 

To construct confidence intervals for 9q, we need to estimate K 4 . Define K n = y l Yrft=i'nt- By 
a simple calculation, we can show that K n — > K 4 a.s. as n —> 00 . Moreover, by letting X n = 
diag(a n , (rf n — l)a^), we can construct a Wald test statistic 


Wn = n(T 0 n - r) T (rs n r T )- 1 (rg n - r) 

to detect the linear null hypothesis Hq : T0o = t, where F E 7Z sx2 is a constant matrix with rank s and 
r E lZ sxl is a constant vector. At the significance level /? E (0,1), we reject Hq if W n > 
where 4 r ^(-) is the cumulative distribution function of Xd- Otherwise, Hq is not rejected. 
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To end this section, we offer some discussions on the model checking. In the context of non¬ 
stationary time series, Ling et al.( 2013) considered two portmanteau tests for non-stationary ARMA 
models, where one is based on the residual autocorrelation functions (ACFs) as in Ljung and Box 
(1978), and the other is based on the squared residual ACFs as in McLeod and Li (1983). However, 


their methods are hard to be implemented for model (1.2). To see it clearly, define the lag -k residual 
ACF of {rjt} as 


Pk ~ 


T!t=k+i(Vt ~ Vn)(Vt-k ~ V* n ) 


where ry* = ^ Ylt=i Vt- By using the fact that 


* = 1 , 2 ,..., 


~ 4>0 — 4>n . , x . /«o" i 

Vt = - 7 ^^sign(?/ t _i) + rjt\ / — and 

VV 


lim — 


n—>oc n 


Vt 


1 a.s., 


t =l 


a simple calculation entails that 


1 


y/n(<t> n — <t>o) 1 

- -j= -V Vt-kSign{y t -i) +o p (l). 

y/ao -n ^ 


Vnpk = —j= Y] VtVt-k 

Wn zJ 

t=k -\-1 v t=k-\-l 

Since { 7 yi_fcsign(yt_i)} is neither a stationary nor martingale difference sequence, it is hard to deter¬ 
mine the limit of r ’ } Y^t=k +1 %-fcSign(yy£_i), and hence that of y/np£. Thus, the classical portmanteau 
tests are not feasible for model and how to check the adequacy of model ( | 1 . 2 [ ) is still a chal¬ 

lenging open question. 


5 Simulation studies 


In this section, we carry out simulation studies to assess the performance of the estimator of 70 , the 
test of stability, and the QMLE of 9q in finite samples. We generate 1000 replications of sample size 
n = 100 and 200 from the following DARWIN(l) model: 


Vt = 0.5yt_i + r]t\J OiQyl_ v 


-3 


where pt is taken as 1V(0,1), the standardized Student’s t$ (sts) with density f(x) = 3 _^ (l+x 2 /3) 
and the Laplace distribution with density f(x) = ^ exp(— y/2\x\), respectively. Here, we set 4>o = 0.5 
and let the value of ao vary corresponding to the cases of 70 > 0 , 70 = 0 , and 70 < 0 , respectively; 
see Tabled) 









Table 1: The values of 70 and a ;? when 0 q = 0.5 is hxed and 07 varies. 



N( 0,1) 



st 5 



Laplace 


ao 

7o 

CT 70 

a 0 

7o 

°70 

a 0 

7o 

4 > 

3.1 

-0.0297 

1.2326 

4.1 

-0.0289 

1.3355 

5.0 

-0.0143 

1.4357 

3.3058 

0.0000 

1.2328 

4.3697 

0.0000 

1.3368 

5.1726 

0.0000 

1.4396 

3.5 

0.0265 

1.2326 

4.5 

0.0133 

1.3374 

5.4 

0.0182 

1.4443 


Table [ 2 ] shows the empirical mean (EM), empirical standard deviation (ESD), and asymptotic 
standard deviation (ASD) of the QMLE 6 n and the Lyapunov exponent estimator 7 ^. The ASD of 
7 n and Q n is calculated by the asymptotic covariance matrix in Theorems 3.1 and 4.1, respectively, 
where the theoretical values of a^ g are given in Table |TJ and the theoretical values of K 4 are 3, 9 
and 6 for A"(0,1), sts, and Laplace distributions, respectively. From this table, we can see that the 
larger the sample size, the closer the EMs and their corresponding true values, and also the closer 
the ESDs and ASDs. Particularly, 7 n performs well even though 70 is very small. To assess the 
overall performance of 7 n , Figj2] plots the empirical density of >/n( j n — 70) for different ao in Table 
[l] From Figj2j we find that 7 n has a good performance in all cases. 


Next, we examine the performance of test statistic T n in (3.1) for testing the hypothesis Hq : 


70 = 0 against H\ : 70 / 0. Fig. [3] shows the power and size of T n for n = 100 and 200 when ao 
varies, and the size of T n corresponds to the case that ao = 3.3058 (for 1V(0,1)), 4.3697 (for sts) and 
5.1726 (for Laplace); see Table [IJ From this figure, we can see that T n has a very precise size, and 
overall, the power to detect the instability is significant, even when sample size is small. 

6 An empirical example 

This section applies the DARWIN(l) model to study the daily exchange rates of New Taiwan Dollars 
(TWD) to United States Dollars (USD) from January 1, 2007 to December 31, 2009, which has in 
total 692 observations. The log-returns of this exchange rate series, denoted by , are plotted 

in Fig|4} 


First, we use model (1.2) with the QMLE estimation to fit {yt} by 


yt = 0.36 66( 0 166 i)yt_i + 7 *-^19.05 86 ( 5 . 3083 ) 


( 6 . 1 ) 
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a 0 

7o 



n = 100 



n = 200 






OL n 

In 


(y. n 

% 

3.1 

-0.0297 

EM 

0.4995 

3.0672 

-0.0308 

0.4966 

3.0845 

-o.o; 



ESD 

0.1721 

0.4442 

0.1138 

0.1226 

0.3182 

0.0' 



ASD 

0.1761 

0.4384 

0.1110 

0.1245 

0.3100 

0.0' 

3.3058 

0 

EM 

0.4978 

3.2943 

-0.0029 

0.5032 

3.2912 

0.01 



ESD 

0.1767 

0.4795 

0.1155 

0.1253 

0.3291 

o.o; 



ASD 

0.1818 

0.4675 

0.1110 

0.1286 

0.3306 

0.0' 

3.5 

0.0265 

EM 

0.5086 

3.4548 

0.0200 

0.4997 

3.4721 

o.o: 



ESD 

0.1866 

0.4877 

0.1143 

0.1292 

0.3464 

0.0' 



ASD 

0.1871 

0.4950 

0.1110 

0.1323 

0.3500 

0.0' 

4.1 

-0.0289 

EM 

0.5013 

4.1054 

-0.0236 

0.4985 

4.1129 

-o.o: 



ESD 

0.2030 

1.0845 

0.1181 

0.1442 

0.9438 

0.0' 



ASD 

0.2025 

1.1597 

0.1156 

0.1432 

0.8200 

o.o; 

4.3697 

0 

EM 

0.5162 

4.3623 

0.0037 

0.4960 

4.4137 

o.oi 



ESD 

0.2127 

1.1575 

0.1160 

0.1497 

0.9874 

0.0' 



ASD 

0.2090 

1.2359 

0.1156 

0.1478 

0.8739 

o.o; 

4.5 

0.0133 

EM 

0.5066 

4.4341 

0.0172 

0.4906 

4.4976 

0.0 



ESD 

0.2115 

1.1814 

0.1235 

0.1463 

0.8349 

o.o; 



ASD 

0.2121 

1.2728 

0.1156 

0.1500 

0.9000 

o.o; 

5.0 

-0.0143 

EM 

0.5050 

5.0073 

-0.0107 

0.4921 

4.9650 

- 0.0 



ESD 

0.2229 

1.0879 

0.1184 

0.1616 

0.7920 

o.o; 



ASD 

0.2236 

1.1180 

0.1198 

0.1581 

0.7906 

o.o; 

5.1726 

0 

EM 

0.4963 

5.1454 

0.0055 

0.5005 

5.1590 

- 0.01 



ESD 

0.2241 

1.1739 

0.1195 

0.1598 

0.8347 

o.o; 



ASD 

0.2274 

1.1566 

0.1200 

0.1608 

0.8179 

o.o; 

5.4 

0.0182 

EM 

0.4962 

5.3933 

0.0148 

0.5023 

5.4004 

0.0 



ESD 

0.2305 

1.2174 

0.1208 

0.1677 

0.8702 

o.o; 



ASD 

0.2324 

1.2075 

0.1202 

0.1643 

0.8538 

o.o; 


m 



rj ~ N( 0 , 1 ) 
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Figure 2: The histograms of o) when r/ is N( 0,1), st 5 , and Laplace distribution, respectively. 

The curves are the densities of 1V(0, <r^ 0 ). Here the true parameter <j>Q = 0.5 and the values of ao and 
cr ^ 0 are given in Table [lj respectively. The top, middel, and bottom panels correspond to 70 < 0, 
70 = 0, and 70 > 0, respectively. The sample size is 200. 
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(a) r| ~ N(0, 1) (b) r| ~ sis (c) r| ~ Laplace 





Figure 3: The power and size of T n at the significance level 5% when rj is N( 0,1), sts and Laplace 
distribution, respectively. 



600 690 


Figure 4: The log-return of daily exchange rates of New Taiwan Dollars (TWD) to United States 
Dollars (USD) from January 1, 2007 to December 31, 2009. 
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where the values in parentheses are estimated standard errors. Based on the residuals { rjt }, Figj5] 


plots the ACF and PACF of {rjt} and {rjt}- From this figure, it seems that model (6.1) is adequate. 
Next, we use the Mira test and the Cabilio-Masaro test in R package lawstat to test the symmetry 
of rjt, and find that their p-values are 0.2985 and 0.2280, respectively. Therefore, we accept the 
hypothesis that rjt is symmetric at the significance level 5%. Moreover, we use the Wald test statistic 
W n to detect the hypothesis Hq : cf >o = 0. The p-value of W n is 0.0273, and it turns out that we 
can reject Hq at the significance level 5%. On the other hand, we find that the estimated Lyapunov 
exponent } n = 0.0001 with = 1.4651, and this implies that the value of test statistic T n = 0.0019. 
Clearly, the null hypothesis of 70 = 0 is not rejected at the significance level 5%. Thus, there is no 
statistical evidence against the hypothesis that the return process is stable, i.e., log volatility of yt is 
a random walk. 
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Figure 5: The top (or bottom) panel is the ACF and PACF of {rjt} (or {rjt}) 


It is also interesting to apply model (1) to fit {yt}- The fitted model with the QMLE estimation 
is 

yt = 0.0317( 0 .o32i)yt-i + 74 ^ 9.3 x lO^-^gQ) + 0 . 3349 ( 0 .o 2 i 4 )l/i-i> (6-2) 

where the values in parentheses are estimated standard errors. Clearly, the estimate of ujq is very 
close to zero. Let {rjt} be the residuals of model (|6.2[). Figj 6 ] plots the ACF and PACF of {rjt} and 
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Figure 6: The top (or bottom) panel is the ACF and PACF of {^} (or {rff}) 


{rtf}, suggesting that model (6.2) is not adequate. Based on these facts, it implies that model (6.1) 


is preferred to fit {yt}- To gain more insight, we plot the estimated log volatilities of models (6.1) 
and (6.2) in Fig jTj The log volatility of model (6.2) is bounded from below by the logarithm of the 
intercept term, while there is no such bound in model (6.1). Among years 2007-2009, the financial 
crisis happened so that the log volatilities probably tend to have no lower bound, and hence this 


might lead to the preference of model (6.1) in fitting {yt}- 


7 Concluding remarks 

This paper proposes a DARWIN(l) model. This new model is non-stationary and heteroskedastic, 
but unlike non-stationary DAR model, it includes a close-to-unit root type behavior of volatility in 
the stability case. In order to study the stability of the DARWIN(l) model, an easy-to-implement 
estimator of the Lyapunov exponent and a test for the stability have been constructed. Moreover, 
this paper studies the QMLE of the DARWIN(l) model, and finds that the asymptotic theory of the 
QMLE is invariant regardless of the stability of the model. An empirical example on TWD/USD 
exchange rates illustrates the importance of the DARWIN(l) model. 

As one natural extension work, one can apply the principle of setting the intercept to zero to 


14 






























0 


Loft volatility in DARW1N(1) model 


-3 



1 100 200 300 400 500 600 690 

l 


-7 

-8 

-9 


-10 


-11 


-12 


Log volatility in DAR(l) model 



1 100 200 300 400 500 600 690 

i 


Figure 7: The top (or bottom) is the log volatility in model (6.1) (or (6.2)). 


high-order DAR models or conditional heteroscedastic models, which would produce a branch of 
models that totally differ from the classical counterparts in the literature. Finally, it is worthy 
noting that how to test the DARWIN model in the null against the counterpart DAR model in the 
alternative is an interesting open question. The traditional tests, e.g., the likelihood ratio test, the 
Lagrange multiplier test and the Wald test, or their one-sided counterparts, may not work due to 
the nonstationarity of DARWIN models. How to develop powerful tests for such hypothesis would 
be a promising but challenging direction for future study. 


Appendix: Technical Proofs 


Proof of Theorem 3.1 


By the definition of 7 n , we have 


1 n _ I I 1 n 

In = -Y'log = - Vlog|c/>o + ?ftsign(y t _i)vao| 

v lyi-iD n “ 1 
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Since (r/isign(yo),r?tsign(y t _i)) = (771 ,..., 77 *) by the induction over t > 1, we have 

d 1 n 

In = - V log Uo + r)tV®o\- 
7=1 

Thus, the result holds. 

Proof of Theorem 4.1 

A simple calculation yields that 

,- n n 

4>n = <t >0 + —— V 77tsign(y t _i) and S n = — V rfi - (<j> n - </> 0 ) 2 . 

n ' n ' 

7=1 7=1 

Note that {? 7 isign(yt_i)} is a martingale difference with respect to Ft = a(m,i < t ). Thus, <f> n is 
unbiased. Also, E(a n ) = ^^ao. Thus, (i) holds. 

By Theorem 2.19 in Hall and Hedye (1980), we have ^ XVi ?? 7 sign(y t _i) ~1 0 a -S., which implies 
that (f> n —> (f >0 a.s., and then a n —> ao a.s. by the strong law of large numbers. Thus, (ii) holds. 

For (iii), note that 

Vn($n ~<t>o) = E™= 1 ? ?7sign(y t _i), 

Vn(a n - a 0 ) = J2t=i(Vt ~ !) “ 7 = (vVV - 0o)) • 

By the Cramer-Wold device and the martingale central limit theorem in Brown (1971), the proof of 
(iii) is trivial. 
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